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result from the second method has a signal significance of 5.4a and is the first observation of 
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Measurements involving heavy vector boson pairs 
(WW, WZ, and ZZ) are important tests of the elec- 
troweak sector of the standard model (SM). Deviations 
of the production cross section from predictions could 
arise from anomalous triple gauge boson interactions [l[ 
or from new resonances decaying to vector bosons. Fur- 
thermore, the topology of diboson events is similar to 
that of events in which a Higgs boson is produced in asso- 
ciation with a W or a Z , allowing diboson measurements 
to provide an important step towards future measure- 
ments of Higgs boson production. 

Diboson production has been observed at the Tevatron 
in channels in which both bosons decay leptonically . 
Extraction of the diboson signal in hadronic channels 
is more challenging because of significantly larger back- 
grounds. In addition, due to limited detector resolu- 
tion, it is difficult to distinguish hadronically decaying 
W bosons from Z bosons. We report on two measure- 
ments of the cross section, a(pp — > WW + WZ), with the 
CDF II detector [1] that use different techniques applied 
to the leptonic decay of one W and the hadronic decay 
of the associated W or Z {WW jW Z — > Ivqq, where £ 
represents a high-p-p electron or muon). Our result rep- 
resents the first observation of this signal in the lepton + 
jets channel. Evidence has previously been reported by 
the DO collaboration and the CDF collaboration set 
a limit on its cross section times branching ratio @. In 
addition the CDF collaboration has reported observation 
of WW + WZ + ZZ in a different hadronic channel with 
large missing transverse energy and jets Q ■ 

The first method uses the invariant mass of the two-jet 
system (Mjj) to extract a signal peak from data corre- 
sponding to 3.9 fb _1 of pp collisions at y/s = 1.96 TeV. 
The second method takes advantage of more kinematic 
information in the event by constructing a discriminant 
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based on calculations of the differential cross sections 
of the signal and background processes. This so-called 
matrix-element (ME) method has been employed in a 
search for a low-mass Higgs produced in association with 
a W boson [§[ and in a measurement of single top produc- 
tion Q. It is expected to achieve greater discriminating 
power and here uses data corresponding to an integrated 
luminosity of 2.7 fb _1 . 

Data samples common to both analyses use trigger se- 
lections requiring a central electron (muon) with Et (pt) 
>18 GeV. The ME method utilizes an additional sample 
derived from a trigger requiring two jets and large miss- 
ing transverse energy (Et) fioj . 

Offline we select events with electron (muon) candi- 
dates with Et (pt) > 20 GeV, and with Et, jet, and 
other kinematic requirements chosen differently for the 
two methods. Jets are clustered using a fixed-cone al- 
gorithm with radius AR = \J (Arf) 2 + (Acf>) 2 = 0.4 and 
their energies are corrected for detector effects [ll| . Cos- 
mic ray and photon conversion candidates are identified 
and removed. 

Further event selection requirements are made to re- 
duce backgrounds and the sensitivity to systematic un- 
certainties. In the Mjj method, we require events to 
have $t > 25 GeV, at least two jets with E T > 15 GeV 
and \rj\ < 2.4, and the dijet vector boson candidate to 
have pt > 40 GeV/c. As a result of these selection cri- 
teria, the Mjj distribution for background is smoothly 
falling in the region where the signal is expected to peak. 
The invariant mass of the dijet vector boson candidate, 
Mjj, is evaluated from the two most energetic jets. Addi- 
tional requirements are made to reduce backgrounds and 
improve the Monte Carlo modeling of event kinematics: 
the transverse mass of the lepton and -St system (Mt(W) 
10]) must be greater than 30 GeV/c 2 , and the two most 
energetic jets must be separated by \Arj\ < 2.5. 

In the ME method, we require events to have Et > 
20 GeV and exactly two jets with Et > 25 GeV and 
\rj\ < 2.0. Additional selection criteria to reduce back- 
grounds and achieve good modeling of the quantities 
used in the matrix element calculation include the re- 
jection of events with either an additional jet of Et > 12 
GeV or a second high-py charged lepton. The latter re- 
duces Z+jets, ti, and leptonic diboson backgrounds. For 
events with an electron candidate, there is a significant 
background from production of multiple jets (multi-jet 
in the following) by quantum chromodynamical (QCD) 
processes, where the electron is faked by a hadronic jet. 
The ME method deals with this background by applying 
stringent selection criteria, while the Mjj method assigns 
a systematic uncertainty to the background shape. The 
reduction of the multi-jet QCD background in the ME 
analysis is achieved by raising the Et cut to 40 GeV, 
requiring Mt(W) > 70 GeV/c 2 , and imposing addi- 
tional cuts on the angles between the jets, the lepton, 
and the Et [H- There 

is a less stringent requirement of 
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Mt(W) > 10 GeV/c 2 imposed on muon events to reduce 
the QCD background in that channel. 

After these selections for both methods, the domi- 
nant background to the diboson signal is a W boson 
produced with accompanying jets (W+jets), where the 
W decays leptonically. Smaller but non-negligible back- 
grounds come from QCD multi-jet (where one jet mimics 
a lepton signature), Z+jets, tt, and single top produc- 
tion. QCD multi-jet events are modeled using data with 
loosened lepton selection criteria. Signal and other back- 
ground processes are modeled using event generators and 
a GEANT-based CDF II detector simulation. The diboson 
signals and the tt and single top back grou nds are simu- 
lated using the PYTHIA event generator 14j . The W+jets 
and Z+jets backgrounds are simulated using the tree- 
level event generator ALPGEN 15j, with an interface to 



> 
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PYTHIA providing parton showering and hadronization. 

The normalization of the Z+jets background is based 
on the measured cross section while for tt and single top 
backgrounds the NLO predicted cross section is used [16| ■ 
The efficiencies for the Z+jets, ti, and single top back- 
grounds are estimated from simulation. The normal- 
ization of the QCD background is estimated by fitting 
the $t spectrum in data to the sum of all contributing 
processes, where the QCD and W+jets normalizations 
float in the fit. In the final signal extractions from both 
methods, the multi-jet QCD background is Gaussian con- 
strained to the result of this£r fit and the W+jets back- 
ground is left unconstrained. 

We now describe the methodology and results from 
each technique. In the Mjj method we extract the signal 
fraction from the data by performing a y 2 fit to the di- 
jet invariant mass spectrum separately for electron and 
muon events. Templates of Mjj distributions are con- 
structed with the multi-jet QCD background, the signal 
WW + WZ processes, and the sum of the electroweak 
backgrounds (Z + jets, W + jets, and tt production). 

Figure Q] shows the fit results superimposed on data 
after the electron and muon samples are combined. Also 
shown is the data Mjj distribution after having sub- 
tracted the estimated background, superimposed on the 
signal Monte Carlo normalized to the fit result. Com- 
bining the two x 2 fit results we get a total of 1079 ± 
232(stat) ± 86(syst) WW/WZ -> Ivjj events, of which 
about 60% are muon events and 40% are electron events. 
The observed significance is 4.6er where 4.9cr is expected. 
The resultant WW + W Z production cross section mea- 
surement is aww+wz = 14.4+ 3.1(stat) ± 2.2(syst) pb. 
The sources of systematic uncertainty in this measure- 
ment are discussed together with those from the ME 
method below. 

In the ME method a probability density P(x) that an 
event was produced by a given process is determined us- 
ing the standard model differential cross section for that 
process. For an event with measured quantities x, we 
integrate the appropriate differential cross section da{y) 
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FIG. 1: Dijet invariant mass distribution of reconstructed 
W/Z — > jj candidates compared to the fitted signal and 
background components (a), and for the corresponding back- 
ground subtracted distribution (b). 



TABLE I: Expected and observed event yields after the ME 
method selection in 2.7 fb -1 of data. 



Process 



Predicted event yield 



WW signal 


446 ± 29 


WZ signal 


79 + 6 


W+jets 


10175 ± 305 


Z+jets 


584 ± 88 


QCD multi-jet 


283 ± 113 


tt + single top 


241 ± 29 


Observed 


11812 



over the partonic quantities y convolved with the parton 
distribution functions (PDFs) and a transfer function: 



P{x) = ^ J da(y)dq 1 dq 2 f(q 1 )f(q 2 )W(y J x). 



(1) 



The PDFs (f(qi) and /(©)) are evaluated according to 
the CTEQ5L parameterization [ljj ■ The transfer func- 
tion W(x, y) relates x to y, encoding the effects of the 
detector resolution. The momenta of electrons, muons, 
and the angles of jets are assumed to be measured exactly 
and a mapping of measured jet energy to partonic energy 
is derived using the full detector simulation. The inte- 
gration is performed over the energy of the partons and 
the longitudinal momentum of the neutrino. The ma- 
trix element is calculated with tree-level diagrams from 
Event probability densities are calcu- 
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lated for the signal processes as well as for W+jets and 
single top background processes. The event probabili- 
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ties are combined into an event probability discriminant: 

EPD — P signal j (P signal "i~ Pbackground) ) where P signal — 
P\VW + PwZ an d Pbackground = Pw+jets + Psingle top- 

We make templates of the EPD for all signal and back- 
ground processes and ultimately extract the signal using 
a fit of the observed EPD distribution to a sum of the 
templates. The expected event yields are as shown in 
Table U for the ME method's event selection. 

Figure [2] shows the dijet mass in bins of EPD. Most 
of the background events have low EPD. Events with 
EPD > 0.25 have a dijet mass peak close to the expected 
W/Z resonance, and the signal-to- background ratio im- 
proves with increasing EPD. 



ure [3] We use a binned-likelihood fit of the observed 
EPD to a sum of templates, testing both a background- 
only hypothesis and a signal-plus-background (s + b) hy- 
pothesis. Systematic uncertainties, discussed further be- 
low, are included in the fit as constrained parameters. 
We perform pseudo experiments to calculate the prob- 
ability (p-value) that the background-only discriminant 
fluctuates up to the observed result (observed p-value) 
and up to the median expected s + b result (expected 
p-value). We observe a p-value of 2.1 x 10~ 7 , corre- 
sponding to a signal significance of 5.4er, where 5.1er 
is expected. The observed WW + WZ cross section is 
aww+wz = 17.7 ± 3.1(stat) ± 2.4(syst) pb. 
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FIG. 3: Observed EPD distribution superimposed on distri- 
bution expected from simulated processes. 



FIG. 2: Mjj for events with (a) EPD < 0.25, (b) 0.25 < 
EPD < 0.5, (c) 0.5 < EPD < 0.75, and (d) EPD > 0.75. 

Before comparing the observed EPD to the prediction, 
we validate the Monte Carlo modeling of the quantities 
that enter the matrix element calculation. We compare 
the observed distributions to the predicted ones in con- 
trol regions with very little signal and also in the signal- 
rich region. The different regions are chosen according 
to the invariant mass of the two-jet system (Mjj): the 
signal-rich region has 55 < Mjj < 120 GeV and the con- 
trol regions cover the rest of the Mjj range. We also 
check the modeling of the properties (mass, px, and 77) 
of the leptonic W boson and the hadronic W or Z boson 
candidate. All of these quantities are well described by 
the simulation for our event selection. There is a small 
discrepancy in the description of Mjj in the control re- 
gions, as is visible in the \aw-EPD region of Figure [2] 
Associated with this discrepancy we assign a systematic 
mismodeling uncertainty which is derived in the control 
regions and extrapolated through the signal region. This 
uncertainty has a negligible effect on the results, because 
most background events lie in the first few bins of the 
EPD distribution. Small changes in modeling of those 
background events do not change the shape of the EPD. 

The observed and predicted EPDs are shown in Fig- 



We consider several sources of systematic uncertainty 
in both methods, taking into account their effect on both 
the signal acceptance and the shape of the background 
and signal templates. The uncertainty on the normaliza- 
tion of the backgrounds is taken as part of the statistical 
uncertainty. In the Mjj method the largest systematic 
uncertainties are due to the modeling of the electroweak 
and QCD shapes, about 8% and 6% respectively. In the 
ME method the uncertainty in the jet energy scale is 
the largest systematic uncertainty, at about 10%, which 
includes contributions both from the signal acceptance 
and from the shapes of the signal templates. In the Mjj 
method this uncertainty is about 6%. Both methods in- 
clude an uncertainty of about 5% due to initial and final 
state radiation and a 6% uncertainty on the integrated 
luminosity. Smaller contributions arise from PDFs, jet 
energy resolution, the factorization and renormalization 
scales used in the VF+jets simulation, and trigger and 
lepton identification efficiencies. 

One measure of how the two methods are correlated 
is the expected overlap of WW + WZ signal. Account- 
ing for the different integrated luminosities used, 15% 
of the signal in the Mjj analysis is common to that in 
the EPD analysis. Conversely, 29% of the signal in 
the EPD analysis is common to that in the Mjj analy- 
sis. This corresponds to a statistical correlation of about 
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21%. If we assume the systematic uncertainties are 100% 
correlated, then the total correlation between the two 
analyses is 49%, leading to a combined [l|| result of 
vww+wz = 16.0±3.3(stat + syst) pb. Because the total 
uncertainties on the two input measurements are similar, 
the combined central value does not depend significantly 
on the correlation assumed. The total uncertainty in the 
combined result increases with increasing correlation and 
we quote the value assuming maximum possible correla- 
tion. The signal overlap with the CDF WW + WZ + ZZ 
observation in the ^T+jets channel [7J is also studied. 
While that analysis requires much larger $t, it does not 
veto events with identified leptons. We found that about 
15% of the WW + WZ signal from the ^r+jets analysis 
appears in the analyses presented here. 

In summary, we observe WW + WZ production in the 
lepton plus jets plus $t final state. We perform two 
searches: one seeking a resonance on top of a smoothly 
falling dijet mass distribution, and another building a dis- 
criminant using a matrix element technique. The com- 
bined WW+WZ cross section from these two methods is 
measured to be <Jww+wz = 16.0±3.3(stat + syst) pb, in 
good agreement with the prediction of 16.1 ± 0.9 pb [2fjj. 
Measurements of these diboson processes are tests of elec- 
troweak theory and a necessary step toward validating 
Higgs boson search techniques at the Tevatron. 
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